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We consider the problem of obtaining the scattering length for a fermion colliding with a dimer, 

formed from a fermion identical to the incident one and another different fermion. This is done in 

, the universal regime where the range of interactions is short enough so that the scattering length a 

I for non identical fermions is the only relevant quantity. This is the generalization to fermions with 

f— ^ different masses of the problem solved long ago by Skorniakov and Ter-Martirosian for particles 

pvi with equal masses. We solve this problem analytically in the two limiting cases where the mass 

of the solitary fermion is very large or very small compared to the mass of the two other identical 

f^ fermions. This is done both for the value of the scattering length and for the function entering the 

Skorniakov- Ter-Martirosian integral equation, for which simple explicit expressions are obtained. 
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I. INTRODUCTION 



C/3 

Cd Recent years have seen very important developments in the field of ultracold fermionic gases [ill. In addition to its 

own intrinsic interest, this field has a strong overlap with problems arising not only in condensed matter physics, where 
the degenerate electron gas is an essential ingredient of metal physics, but also in nuclear physics and in quark matter 
physics [2]. A major strong point of these cold gases systems is the extreme simplicity of the effective interaction 
r^ between fermions of different species, which can be fully described by the mere knowledge of the scattering length a, 
^H while there is no effective interaction between same species fermions, since this has to occur in the s-wave channel 
(_J because of the very low temperature and energies involved and this s-wave scattering is forbidden by Pauli principle 
Ci (we assume that no p-wave or higher- wave resonance occurs). In addition the scattering length can be changed, in 

H the vicinity of Feshbach resonances, by merely modifying the static uniform magnetic field applied to the trapped 
I gas. This offers the fascinating possibility to control at will the strength of the interaction. 

The first explorations of these gases have dealt with systems containing two different kinds of fermions with same 
particle number and same mass, in practice fermionic atoms belonging to the same element but being in different 
O hyperfine states. This is the most convenient and the most natural system to explore experimentally, and in addition 
the physical situation is the same as for electron in metals which may have their spins up or down and are roughly 
in equal number. A fair exploration has been made of the normal state and mostly of the superfluid state where the 
whole BEC-BCS crossover has been studied to a large extent. More recently the research activity has moved toward 
1/-V systems where there is a sizeable imbalance between the two fermionic populations ^Ij and where a great deal of 

Q^ interesting physics has been found. Quite recently some attention for fermions with different masses has arisen. This 
00 is because new interesting physics, both for normal and superfluid phases, is expected when the mass ratio becomes 
C^ large [3]-(n] • Experimental developments in this direction are in progress in several laboratories in the world (TUl [TT] 
r — and degeneracy has already been reached in these mixtures. Similarly several quite recent theoretical works have 
f^ focused on this situation [121 I13j , with an emphasis on systems with a small number of particles where multiparticle 
^^ bound states are expected to occur. 

^^ As we have mentioned the most interesting situations occur when there is a s-wave Feshbach resonance, which is 

^ directly linked to the appearance of a bound state, or molecular state, or dimer, between two atoms belonging to 

different fermion species. This will occur on the so-called BEC side of the resonance, and on this side these dimers 
play clearly an essential role in the physics of the gas. As it is quite common we call the two involved fermion species 
^ t and ^ spin, by analogy with electrons in solids, even if we are actually interested in fermions belonging to different 
atomic elements. If the masses of these atoms are respectively m^ and to^, the mass of the dimer is the total mass 
M = m^ + mi. The binding energy is Ej, = l/(2/.ta^) where /i = m^mi/{m^ + mi) is the reduced mass. In addition 
to these dimers we have naturally in general on the BEC side isolated m-^ and mi fermionic atoms. This will occur 
for example at zero temperature if there is an imbalance between the two fermionic populations so that all fermions 
can not be paired into dimers. Or this may also occur at non zero temperature because the dimers will be thermally 
broken. Hence scattering between dimer and isolated atoms will happen, and in cold gases this scattering will be 
fully characterized by the atom-dimer scattering length 03. Clearly this is, beyond the up-down scattering length a, 
one of the basic quantities coming in the physics of these systems. This is the quantity we are interested in. Quite 
specifically we will investigate the scattering length 03 of a f atom by a dimer. 



Naturally an important feature in these kind of systems is their stability [T]. It is known that inelastic three-body 
collisions are a major process for loss rate in these systems, being responsible for their limited experimental lifetime. 
Obviously the existence of dimers favours these processes since two atoms are already fairly nearby in these bound 
states. This is indeed known to be a major problem in bosonic systems. Fortunately in fermionic systems, with which 
we are concerned here, these processes are comparatively restricted by the fact that, qualitatively, Pauli exclusion keeps 
the same species fermions apart, limiting the overlap with deep bound states responsible for relaxation. This is origin 
of the fairly long lifetime found for fermionic systems made of the same atomic element. However the situation might 
be somewhat different for fermions with different masses and this is obviously a question under current experimental 
investigation. Moreover in these systems dimer-dimer collisions provide also an important relaxation channel [T]. 
Finally [13] three-body and four-body resonances occuring in these systems might lead to efficient additional loss 
channels. In the present paper we will assume that, nevertheless, the stability of the system is high enough for our 
investigation to be meaningful. 

The calculation of 03 is a long standing problem and, for equal masses, the problem has been solved long ago by 
Skorniakov and Ter-Martirosian [14]. The generalization to different masses is naturally quite easy, either following 
the original paper [nj or making use of our recent work [IP , as it has been done recently by Iskin and Sa de Melo 
[16j . The resulting equation has been quite recently fully investigated numerically by Iskin 17 . Another approach 
has been used by Petrov |18j who has obtained 03 by solving the 3-body Schrodinger equation with the appropriate 
boundary conditions imposed by the scattering length a. 

While the above works solve completely the problem as far as numbers are concerned, we feel that a deeper 
understanding of the involved equations is warranted for such a basic problem. Moreover such a full control provides 
a first step to address more complicated problems which will arise in many-body problems, involving the same building 
blocks, which can be formulated with the same Green's functions formalism, using the same Feynman diagrams as we 
have done for example in [19) . Such a control means in principle a complete analytical investigation. Unfortunately 
this does not seem possible in the general case. On the other hand we have been able to carry out this investigation 
in the limiting cases where the mass ratio r = m^jm^ is either very large or very small and the purpose of this paper 
is to present our results. We note that for the small r limit Petrov [TS^ has done briefly a similar investigation based 
on a Born-Oppenheimer approximation within the Schrodinger equation. Here we address this problem within the 
Skorniakov and Ter-Martirosian approach, in a more systematic way and with more specific and precise results. After 
presenting the basic equation in section |ll| we consider in section |III| the case where the mass ratio r is very large, 
and then in section [TV] the case where it is very small. Since this last case is more involved section [TV] contains several 
subsections. Let us finally mention that an important ingredient in all the few-body problems is the existence of 
bound states of the whole system under study, in our case bound states of the t atom and the dimer. This has been 
studied in great details quite recently by Kartavtsev and Malykh [20] and these bound states include in particular the 
Efimov states [21] which have received recently a great deal of attention. However since we are here only interested 
in fermions and in s-wave states, these bound states do not come in our problem since they appear only for non-zero 
angular momentum [20] . 



II. BASIC EQUATIONS 

Our basic equation is obtained by writing directly [15^ the integral equation for the full dimer-fermion scattering 
vertex T3{pi,P2tP), where P is the (conserved) momentum-energy of the particles P = {P,i?}, with P the total 
momentum and E the total frequency of incoming particles, pi is the momentum-energy of the incoming ^ fcrmion 
and p2 its outgoing value. 

The scattering length is obtained by considering the case of total momentum P = and total energy E = —Et,. 
The incoming and outgoing momentum-energy of the scattering fermion should also be zero pi = P2 = 0. Introducing 
the "on-the-shell" value of T3 for the incoming fermion and the related function a3(p) according to: 
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the atom-dimer scattering length 03 is merely given by 03 — 03(0). Here fix is the atom-dimer reduced mass 
^T = m^M/{m^ + M) = m^{rn^ -f ?Ti^)/(2m^ + m^). Then the integral equation for T3, restricted to "on-the-shell" 
values, becomes the generalization for different masses of the equation of Skorniakov and Ter-Martirosian for 03 (p), 
namely: 
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where we have used reduced units by setting a-^/a = 03 and p = p/a and q = q/a. The angular integration is easily 
performed, since 03(5) depends clearly only on |p|. Using the simpler notations R = ^jl/^it^ R' = 2/i/m^ (they are 
related by i? + i?'^/4 = 1) , together with x — \p\ — p, y — \c^\ — q and f{x) — 03(5), one finds: 

1 + Vi + Rx^ 1 + x^ ^R' Jo xy 1 + a;2 + j;2 _ jif^y y ' 

For equal masses, R — 3/4 and R' = 1 and solving this equation gives the well-known result [T3] for the atom-dimer 
scattering length 03 = 1.18 a. 

The generalization of the Skorniakov and Ter-Martirosian equation following the above lines has been obtained by 
Iskin and Sa de Melo [T5], who have solved the equation numerically and plotted a^/a as a function of the mass ratio 
r = mi/m^. This has been taken up in great details by Iskin |17) . Let us recall that the homogeneous equation 
corresponding to Eq.(l3]) has no solution, corresponding to the physical fact that [20l [21] there are no Efimov states, 
nor any other bound state, for the three fermions problem we are considering. 

Before turning to the limiting cases we are interested in, we present first some general results which will be useful in 
the following. Let us first consider the behaviour of f{x) when a; — >■ cxi. If one assumes (which is confirmed numerically) 
that, in the integral in the right-hand side of Eq.dsl), f{y) goes rapidly enough to zero for large y (possibly with an 
oscillating behaviour) so that the integral converges, one can consider y as effectively bounded, and for x — ^ 00 expand 
the logarithm. This leads from Eq.(|3|, in this large x range, to the simplified equation: 
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Apparently this leads to f{x) ^ 1/x, in contradiction with our starting hypothesis since, in this case, the integral in 
right-hand side of Eq.Q diverges. The only possible escape is that the coefficient of 1/x^ in this right-hand side is 
zero. This leads to the simple relation: 

dy m = 2 (^) 

This relation is well satisfied numerically for all the mass ratio r that we have investigated. 
Another simple relation is obtained by merely setting a; = in Eq.([3]). This gives: 

2 TT Jo l + y^ 

III. VERY HEAVY MASS 

In the case where the down-spin particle has a very heavy mass, one has a situation quite close to the case where 
this mass is infinite m^ = 00. In this last case, the physical situation is simpler than in the general case. Indeed the 
down particle acts merely as a fixed diffusion center for the two up-spin particles. Numerical calculations give a^ = a 
with excellent precision in this case. 

Indeed analytically, since we have then R = I and R' -^ 0, Eq.Q reduces in this case to: 

fix) 1 2 p 1 .„. 
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and we have found that it has the simple, but non-trivial, solution: 
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Hence we find a^ — f{0) — I, which shows that the above conclusion from the numerical results is actually exact. We 
can check on this simple solution that Eq.([5]) (and naturally also Eq.Q) is indeed satisfied. 



This simple but non trivial result can be easily understood if we go back to the r space formulation and consider 
the wave function of the system. In the limit of a very heavy mass ?7^^ = oo, the spin down particle does not move and 
it acts as a simple diffusion center, located at the origin for convenience, for the two spin up particles. Since they do 
not interact, the eigenstates are properly antisymmetrized products of single spin up particle wave functions. When 
one spin up particle is far away from the origin, we want to describe a physical situation where one spin down particle 
is in the bound state linked to the diffusion center, described by the (non normalized) wave function exp(— ri/a)/ri. 
The other spin up is in a scattering state corresponding to an incident particle described by a plane wave exp(iko.r2) 
and a scattered particle in a diverging wave exp(iko.r2)/r2. We are only interested by the s-wave component of this 
wave function, which is sin(fcor2 + 5o{ko)) / {kor2) , with So{ko) the s-wave phase shift. This is the boundary condition 
for our scattering problem, but since the spin up particles do not interact this wave function is an exact eigenstate 
for any value of r2. The full antisymmetric wave function is then: 

_ exp(-ri/a) sin(fcor2 +5o(fco)) exp(-r2/a) sin(fcori + (5o(fco)) ,_,^. 

ri kor2 r2 k^ri 

It is already clear at this stage that the scattering amplitude is the same as the one of a single spin up particle scattered 
by the diffusion center. The spin up particle in the bound state does not play any role since the two particles do not 
interact. 

If we consider the fcp — > limit, which is the one of interest for us for the scattering length, taking (5o(fco) — —k^a 
into account, we have merely 

vl.(ri,r2) = ^'^PhllM (l-^\- ^M-r2la) L _ a \ ^^^^ 

ri \ r2j r2 V '"i / 

Following Skorniakov and Ter-Matirossian [14], we define a scattering wave function x{t) = liniri^o (''i^lri, r)). 



From Eq.(ll|, we find 



X(r) = 1 - " + °; exp(-r/a) (12) 
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The large distance behavior of x{'>') determines the atom-dimer scattering length as according to 

X(r) « 1 - y (13) 



valid for r ^ a. Looking at Eq.( 12 ) in the limit r ^ a, we find indeed as = a. 

Let us define x(p) as the Fourier transform of xi''')- We know [M] that a3(p) and x(p) are related to each other 
according to 

X(p) = (27r)3 5(p)-4^"^ (14) 

P 

If we Fourier transform Eq.( 12 1, we find x(p) — (27r)'^ (5(p) — 47r ^. _2! 2ip2 - This corresponds to 



"3(P) = ( -2 I ^2^ (15) 
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which is precisely what we found (translated in non reduced units) as a solution of our integral equation. 
The normalisation condition Eq.(l5]) can also easily be understood in this limit. Using Eq.(14), we easily find 



^^"^(^) ^ |-f/(^^(P) 



^-^X{r = 0) (16) 



From Eq.(12| we can directly check that %('' = 0) = 0, and therefore /„ dpa^{p) — 7r/2. However this result is 
obvious physically and therefore quite general. It states that, when the two up spin particles are both at the origin 
(ri = r2 = 0), the wave function 'I'(ri,r2) is zero due to Pauli exclusion principle. Hence it is quite likely that the 
same explanation for Eq.(|5|) holds in the general case where m^ is not infinite, and we show explicitly in Appendix [Al 



that this is indeed the case. 



IV. VERY LIGHT MASS 

This situation turns out to be much more comphcated than the preceding one, and it is the essential focus of our 
paper. 

A. Case mi — 

Let us first investigate the hmit m^ — > 0. We will find a fairly singular result, but this will allow us to understand 
how to handle the case of a very light mass. For m^ -^ 0, we have i? — > and R' -^ 2, so that Eq.(l3]) becomes: 

Quite generally, from Eq.®, f{x) is an even function of x. This allows us to rewrite Eq.(17l as: 

dyMin[i + (:,-y)2]^2^--^ (18) 

; y 1 + X^ 

This convolution product can more easily be handled by integrating by parts. Introducing: 

Hv) ^ r dy^ (19) 

J-oc y 

which satisfies J-"(+oo) = by parity, and is even with respect to y, by parts integration gives: 

/oo 
_^d.^(x-.)^ = -.^ (20) 

This is easily solved by Fourier transform. The Fourier transfom of x/(l + x"^) is: 

t(fc) = -i7rsign(fc) exp(-|fc|) (21) 

but since it appears on both sides of the equation, the result for the Fourier transform J-{k) (which is real and even): 



F{k) = / dyF{y) exp(-iA:y) (22) 

^— oo 

is immediately seen to be: 

-F(fc) = -7r (23) 

Going back to our original function this translates into: 

f{x):^~-Kx5'{x) (24) 

where the derivative 5' of the Dirac 5 function is defined as usual by /_ dx5'{x) F{x) — — F'(0). Indeed inserting 
the above result for f{x) into Eq.(18), it is immediately seen that it is solution. We can check that relation Eq.(l5]) is 
satisfied since: 
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Actually this solution is not unique since, from Eq.(17), one may add any function extremely localized around the 



origin (i.e. a (5-like function), and having a zero total weight. An example of such a function is h(x) + x8' (x). As a 



result ■Kb{x) is also solution, as it is easily directly checked in Eq.(17). This singular degeneracy will disappear when 
we consider in the next subsection a very small, but non zero, rn^ and we look for a regular solution. It makes also 
pointless to discuss the value of the scattering length, although the fact that -Kbix) is solution makes reasonable the 
fact that it is infinite in the present limit. 

If we take Eq.(|6| we see that the right-hand side is zero. This looks to be coherent with the fact that i? = in the 
case we consider. But since we suspect that /(O) goes to infinity, this shows more precisely that the scattering length 
does not diverge as fast as 1/i? ~ m^/(2TOjJ. This will be confirmed in the next subsection. 



Finally we note that the solution Eq.(24) we have found is completely concentrated at cc = 0. Accordingly we 
expect that, for the case of very small m^ which we investigate in the next subsection, the solution will be strongly 
concentrated around a; = and decay very rapidly for large x. We remark as a consequence that the condition we 
required for Eq.(l5]) to be valid is satisfied. 



Very small mi 



Just as in the preceding case we will use a Fourier transform, so it is useful to translate the various quantities and 
relations we are interested in. First of all we have for the scattering length: 
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(to avoid divergences, fc-integrals must be regularized by putting a factor exp(— r/|fc|) and then letting rj — > 0-|_) so 
that: 
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since ^{k) is an even function of k. In the case of the preceding subsection with m^ = 0, this gives indeed 03 = oo 
because we had J'{k) — —tt. 

Similarly the relation Eq.(l5]) becomes: 
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which gives: 
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Let us now expand Eq.(|3]) to lowest order in r = mi/m^, with i? ~ 2r and R' ~ 2(1 — r). The left hand-side 
becomes R f{x)/2. One may wonder if neglecting entirely the term Rx^ is a proper approximation, since even when 
R is small, it may become very large if x goes to infinity. However we have stressed that, for small m^, f{x) is very 
strongly peaked around x = Q and becomes very small for large x. Hence what happens for large x is irrelevant and 
our approximation is legitimate. In the right-hand side we expand the logarithm to first order in i?' — 2 = — 2r. Just 
as in subsection IV A we take advantage of the even parity of f{x) to extend the integral from —00 to 00. After 
multiplication by ttR'x, we obtain: 
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Again in the last term the convolution product will give a simple product in Fourier transform. The Fourier 
transform of l/(l-|-a;^) is merely isign(fc) t{k). Taking into account that the Fourier transform of /(a;) is —d{kF{k))/dk 
from Eq.(27l, we obtain: 
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where we have restricted ourselves to the case fc > 0, which is enough since T{k) is even in k. We see immediately 
that, if we assumed T{k) to be a constant, the result would still be J-{k) — — tt, unchanged with respect to the 
case mi = 0. Introducing accordingly G{k) = J-{k) + tt, performing the derivations and making use of the explicit 
expression of i(fc), we find the homogeneous differential equation: 



[fc(l - e^) G" + (2(1 - e'') - fc) G'] + G = 



(33) 



We see that condition Eq.(30l, which becomes G(0) 



the vicinity of fc = 0, the two independent solutions of Eq.(|33|) behave as k 
solution with a = — 1 



0, is automatically satisfied by Eq.(33|. More precisely, in 
with a 



= -1 ± ^1 + 1/r. Only the 



■\/l -I- 1/r satisfies Eq. ( 30 ) , which provides us with a boundary condition. 



Since r is small, this equation is similar to a Schrodinger equation for which the WKB approximation could be 
applied. It is convenient to make a change of function to bring it precisely under the form of a Schrodinger equation. 
This is done by setting: 
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which leads to: 



Just as G{k), g{k) must satisfy: 
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On the other hand, because of the term e'' in the coefficient of g" , as soon as k is of order of a few units, this coefficient 
will be large compared to the coefficient of g. More specifically, when k > ln(l/r), Eq.(35) will reduce to g" — (we 
might worry that, in the coefficient of g, —k/2 will dominate l/r for large k; however when this occurs we will be 
anyway fully in the regime where g" = is valid, so this does not change our conclusion). In this regime the general 
solution is naturally g{k) = Ak + B. However from Eq.(28) we know that, for large k, T{k) ~ —na^/k. This implies 
G{k) ~ TT — -na^/k and from Eq.(34) g{k) ~ Trfc — 7103. Hence we know that, for large fc, g{k) has an asymptote with 
slope A ~ TT. This fixes the coefficient of the solution of the homogeneous equation Eq.(|35]) satisfying g{0) = 0. Hence 
we have a unique solution as it should be, and the constant B — — ttos will gives us the scattering length. Equivalently 
we can say that, if we take a solution satisfying g{0) = 0, its asymptotic behaviour for large k is g(k) ~ Ak + B and 
we have 03 — —B/A. 



Now, with the above discussion, we can further simplify Eq. ( 35 ) taking into account that we are only interested in 
this equation when r is small. In this case, when k is not very large, the coefficient of g is essentially equal to l/r. 
On the other hand when k is large enough so that this is no longer valid, the equation reduces to g" = 0. Hence in 
all cases Eq.(35) may be simplified into: 
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where we have set r 



WKB solution 



Since we are in a situation where e^ is very small, we can solve Eq.(37) by the WKB method. The solution which 
satisfies g{Q) = is, within a multiplicative constant: 



m 



1 



[Q{k)Y" 



exp 



dx v Q(a 



(38) 



where the arbitrary choice of the upper bound of the int egra l modifies only the arbitrary multiplicative constant. We 
note that, for fc — > 0, Q{k) ~ 1/fc^, and the solution Eq.(38l behaves as g{k) ^ fca + i. This is naturally in agreement 
with what can be obtained directly from Eq.(37), and this corresponds to what has been found for Eq.(33) for small 
e, taking the change Eq.(34) into account. 



Unfortunately we can not extract from Eq.(38) the asymptotic behaviour of g{k) for large k because in this limit 
Q{k) — >■ and we have to deal wi th a turning point for our Eq.(37l which is located at fc = 00. Nevertheless for large 
k we have Q{k) ~ e~^ /k and Eq.(37) becomes: 
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and its solution Eq.(38) simplifies somewhat into: 



g{k) = k^/^ e^l^ exp 



dx 
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As for the standard treatment of the Schrodinger equation within the WKB approximation, we have to solve directly 
Eq. ( 39 1 near the turning point and match the solution to Eq. ( 40 ) . Since the turning point is at infinity it is more 



convenient to bring it to zero by the additional change of variable: 
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which transforms Eq.(39) into an equation for g{z) 



— g{k), namely: 
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with L{z) — 2e"' ln(2/z) and the solution Eq.(40l reads now within an arbitrary multiplicative constant: 
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In contrast with the standard treatment of the turning point in the Schrodinger equation, where the solutions are the 
Airy functions with known mathematical properties, we do not know the solutions of Eq.(42| and we would have to 



find their mathematical properties. However we may notice (and this could be a first step in a mathematical study) 
that, for small 2;, L{z) is a quite slowly varying function of z, compared to z or z'^. In a first step, we may treat it as 
a constant L{z) ~ L. Then the further change of variable z = \/Lx brings Eq.(42) into: 
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for g{x) — g{z). The solutions of this equation are [27 the well known Bessel functions Io{x) and Ko{x). The form 
of the solution Eq.(43) for large z, which becomes now within an arbitrary multiplicative constant: 
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serves as a boundary condition to find the p rop er solution of Eq.(44). Since for large x, we have loix) ~ e^/v27ra; 
while Ko{x) ~ -yAT^e"^, we see from Eq.(45l that the solution we are looking for is just g{x) — Ko{x). On the 
other hand it is known [22] that, for small x,Ko{x) — — Ina; + ln2 — C where C — 0.577215... is the Euler constant. 
Going back to the variable k through x = 2e~'^'^/evfc, we obtain: 



g{k) ~ -+lne 



C+-\nk 



(46) 



The last term, proportional to Infc, is not consistent with Eq.( |39[ ) and its presence is due to our above approximation 
of treating the logarithm L(z) as essentially constant. However if we stay coherent with this approximation and 
notice that the domain we consider for the behaviour of Ko{x) is x ^ 1, we have to dominant order in this domain 
k ~ ln(l/e^) and we have to replace coherently Infc by ln(ln(l/e^). As we have seen 03 is obtained from the ratio of 
the constant to the coefficient of k in Eq.(46|. From e — ^/rniJrri^ this leads to: 



as ~ ln(m^/TO;) - ln(ln(rn^/TO;)) + 2C 



(47) 



Hence to dominant order we find 03 ~ ln(m-^/r?T,^), which is in agreement with the estimate by Petrov [18 in the large 
m-f/m^ regime. However our handling of the Infc term as a quasi constant may look somewhat uncertain and one 



may be more uncertain about the last two terms of our result Eq.(47). Therefore we present in the next section an 
alternative derivation of the scattering length from Eq.(37) which will give us more confidence in our result Eq.(47) 
and will even allow us to go beyond this result. 



D. Alternative derivation of the scattering length 



In the preceding section we have obtained our result by matching the WKB solution at infinity, which requires 
a somewhat uneasy handling. Here, in order to obtain the scattering length, we use another method which avoids 
this problem, or rather puts it on a much more secure basis. It makes use of the fact that Q{k) goes very rapidly 
to zero when k goes to infinity. Hence we can integrate Eq.(37), making use of the boundary conditions 17(0) = 
and g'{0) — which results directly from the g{k) ^ A;3+= behaviour we have found in the preceding section. A first 
integration gives: 



eV(fc)= / dxQ{x)g{x) 
Jo 



(48) 



Another integration, followed by a by parts integration leads to: 



g{k) ^ k I dxQ{x)g{x) 



dx X Q{x) g{x) 



(49) 



Since Q{k) goes very rapidly to zero when k goes to infinity, whereas we know that g{k) — Ak + B for large fc, we 
may replace the upper bounds of the integrals by +00 as soon as k is large. This gives in this large k regime: 



e g{k) — k dx Q{x) g{x) dx x Q{x) g{x) — Ak — B 

Jo Jo 



(50) 



which displays precisely the expected asymptotic behaviour. The scattering length is then obtained from 03 ~ B/A. 
At first this seems useless since the integrals imply the knowledge g{k). However, since Q{k) goes very rapidly to 
zero, we can use one of the approximate solutions which have been discussed above and which would be valid for 
the whole range of integration. A first idea is to use the WKB solution Eq.(38) which leads, within an unimportant 
multiplicative constant, to: 



A = 
B = 



dk [Q{k)] 



3/4 



exp 



dx \/Q{x) 



dkk [Q{k)Y' exp 



3/4 



dx \/Q{x) 



(51) 
(52) 



Since the WKB solution increases rapidly with k while Q{k) decreases rapidly, a further idea is to use a saddle point 
estimate. For B for example the saddle point fcg corresponds to a maximum of the integrant of Eq.(52), or rather its 
logarithm. This gives for fcp the equation: 



4V^ 



(53) 



where we have already used the fact that /cq is large. To dominant order this gives fco — 21n(l/e) — \n.{m^/'m\,). To 
dominant order the saddle point for A is located at the same place, because A and B differ only by the factor k in 
the integrants. For the same reason we obtain to dominant order 03 = B/A ~ fcg ~ ln(77i^/r?T,^), which allows already 
to recover the dominant order of our result Eq. ( 47 ) . Unfortunately this result for fco shows that the saddle point is 



located at a fc value where the argument of the exponential in the WKB solution Eq.(38l is of order one. This implies 



that the WKB solution is no longer reliable because this fco is too large and we are already not far enough from the 
turning point located at infinity. 



Fortunately we have found an approximate solution of Eq. ( 37 1 which matches the WKB solution to the solution at 
infinity, namely ^^^0(0;) with x = 1e~^''^ j t\fk. This solution is only approximate because we have treated fc as a quasi 
constant, compared to e~^ . However the fact that the integrants for A and B are sharply peaked, allowing as we have 
seen the use of a saddle point evaluation, makes this treatment of fc as a quasi constant a very good approximation 
for our purpose. In particular the sharp decrease of Q(fc) makes that we do not have to worry about very large values 
offc. 

Performing the change of variable x = 2e^''/^/e\/fc in A and B, as given in Eq.(50l, with fc quasi constant, we find: 



A = — / dxx Kn (x) 
2 Jo 

B — — dxxk{x)Ka{x) 

2 Jo 

where k{x) is the inverse function oi x = 2e^^l'^ je^fk. It satisfies: 

k(x) = In — - 2 In a; + 2 In 2 - In k{x) 

Actually we do not need a saddle point calculation since we have the exact results 

dx xKq{x) = 1 
dx X \q.xKq{x) = ln2 — C 



f 

Jo 



(54) 
(55) 

(56) 

(57) 
(58) 
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Exact solution of Eq. (3) 

Solution of Eq. (37) 

ln(m^/mi) 

Inlm^/mn — ln(ln(?ri-f/m|)) 

ln(mfjmi) — ln(ln(m^/m|)) + 2C 

Solution of Eqs. (60) and (61) ...••■■' 



4 6 

ln(m|/m;) 



10 



FIG. 1: (color online). Reduced fermion-dimer scattering length 03 — a^/a as a function of the mass ratio 7n-|-/mj_. Full line: 
exact solution from Eq.(|3|. Dashed line: exact solution of Eq.(37l corresponding to the first order approximation of Eq.® 
for small mi/irif. Dotted line: dominant result ln(7Tif/rra^). S mall dotted line: inclusion of the subdominant contribution 
-ln(ln(m-f^/r7i^)). Dashed-dotted line: full asymptotic result Eq.(47l. Dashed-double-dotted line: parametric representation 



Eq.(60l and Eq.(61l 



Since in the integral we have x ^ 1 because -^^0(2;) decreases exponentially for large x, we have to dominant o rde r 
k{x) ~ ln(l/e^) and we may replace k{x) by this approximate value in the last term in the right-hand side of Eq.(56|. 
This leads us exactly again to Eq.(|47]) for the expression of 0.3. 

We may even try to improve on this result by refining our handling of the last term in Eq.(56). We could just 
iterate our procedure. However this would produce new integrals to evaluate and there is no justification to go to 
such a refinement. On the other hand, since we have seen that a; ~ 1 for the domina nt r ange of integration, we may 
just set Ina; ~ 0. This leads to replace, in the last term in the right-hand side of Eq.(56), k{x) by ko- This gives for 
ko the equation: 



ln(r7i-|-/?7i^) — ko + In fco — 2 In 2 



(59) 



while as is still obtained from Eq.(56) by replacing k{x) by fco in the last term in the right-hand side, which yields: 



as = ln(TO-|-/TO^) — In fco -f 2C 



(60) 



This gives us a parametric representation of 03 in terms of ln(?7i^/7n^), with fco as parameter. This is easily plotted. 
The result has the pleasant feature of eliminating the divergent behaviour of Eq.(47l for small ln(m-|-/m^) (which is 
nevertheless not at all disturbing since we are exploring actually the large \n{m^/mi) domain). However the result 
is slightly less satisfactory for large \ii(m^/mj^). This can be improved by noticing that replacing In a; by is just a 
first reasonable try, but that the actual result is rather obtained by replacing it by a constant of order 1. Hence if we 
replace Eq.(59) by: 



ln{m^/m^) = fco + In fco 



(61) 



and treat F as a fitting parameter, we may hope to obtain a better result. Indeed we obtain an excellent result over 
the whole range of variation of ln(?Ti-|-/?TT,^) for F = 2.7 as it can be seen in Fig. fl] 

E. Comparison with numerical results 



As we have mentioned at the beginning, it is naturally easy to find at any stage the numerical answers to the 
questions we raise. For example we display in Fig.nl the numerical solution of the exact problem Eq.([3]) giving 0,3 as a 
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k/a 



FIG. 2: (color online). Solution of Eq.(62l for m^/rai — 10^. To be explicit we use non reduced units for k, which has the 
dimension of a length and is expressed naturally in terms of the scattering length a. Note the very flat part for small k/a. 



function of ln(77i-|-/?Ti^). We can also solve numerically for 03 the equation Eq.(37) we have obtamed for large m-^/m^. 
The result is also displayed in Fig.[T] As it could have been expected the difference between the exact solution and this 
approximate solution of large m^jm^ becomes quite small for ln(m^/m^) ~ 2, i.e. m^ /m\^ ~ 7. We have then shown 
the lowest order approximation a^ — ln(r7T,-|-/TO^) which is a straight line on this graph and which is clearly not so 
satisfactory. Next we have plotted the first two terms of Eq.(47|. It is clear that the — ln(ln(m^/m^)) corrective term 



leads to the proper shape for large ra^/m^. Finally plotting me whole Eq.(47) shows that the constant 2C leads to a 
perfect agreement with the exact result for large m^/m^. However for low m^/m^ this result departs strongly both 
from the exact one and from the exact solution of equation Eq.( |37|), d ue to the divergence of the — ln(ln(7Ti^/7TT,^)) 
term. This deficiency is eliminated by the parametric solution Eq.( |60[ ) and Eq.(61), which gives an agreement with 
the solution of Eq.(37) which almost perfect over the whole range rw^m^ > 1. 



F. The function 03 (p) 



Finally, it is of interest to consider the app rox imate solution 03 (p), or f{x) with our simpler notations, we find for 
the basic integral equation Eq.S. From Eq.(|27| its Fourier transform f{k) is: 



m = -:ruiknk)) = :jrA-k- 



dk 



dk 



m 



(62) 



Vl-e- 

which is easily calculated numerically. An example is given in Fig. [2] It is peaked around fc = and has a width of 
order fco ~ \n{m^/m\J) as it can be seen from the WKB expression Eq.(38) for g{k). Hence f{x) is peaked around 
X — Q with a typical width l/ln(m^/m^) and is accordingly very narrow, as expected. One checks that Eq.(62) gives 
/(O) = 03 from the asymptotic behaviour of g{k). The large x is behaviour is obtained from the small k dependence 
of /(fc), which is /(fc) ~ tt — k~^/^g{k) with g{k) ^ fc^ + i as we have seen. This leads to an extremely fast decrease 
fix) ~ l/a;i+i/^ at large x. 

However we are essentially interested in the domain where f{x) is clearly different from zero, which is obtained for 
fairly small val ues of x since f{x) is narrow. These will come from /(fc) with fairly large values of fc. In this range we 



can set in Eq.(62| 1 — e"*^ ~ 1 and take the approximate solution ^(fc) = 27r/i'o(2e~'^/^/(e-\/fc)) where the prefactor 
provid es t he correct asymptotic behaviour for g{k). On the other hand if a; ^ the first term in the right-hand side 
of Eq.(62) gives no contribution. Taking into account the even parity of /(fc), integrating by parts and taking the 



complex conjugate we obtain; 



fix) ^2xlm / dke-'^'^Koi 



ty/k ■ 



(63) 
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FIG. 3: (color online). For m^/mi — 10 corresponding to 6 ~ 3.08, full line: Eq.(64l, dotted line: Eq. (671, dashed-dotted line: 
exact numerical solution of Eq.(31 1. 



For small k the Bessel function Kq goes very rapidly to zero, while for very large k the Besscl function varies slowly 
while the oscillating factor e*''^ leads to destructive interferences. Hence the dominant contribution arises when the 
argument of the Bessel function is of order unity, that is again when k is of order fcg — lii(?7i-|-/?7i^) = 2 ln(l/e). Hence 
we can again treat in the argument of the Bessel function fc as a quasi constant and write ^/k ~ •\/21n(l/e). Taking 
t = 2e"'''/V(e\/21n(l/e)) as new variable, and making use of [22] J^ dtt^' Ko{t) = 2^'-^ r2((l + ^)/2), where r{x) is 
the Euler function, we obtain: 



/(x) = X Im I [{mi/m^) ln(TO^/m;)]*^ T^ (ix) | 



(64) 



Since [22], for z — )• 0, we have T{z) ~ 1/z 
coherence of our approximate treatment. 



C, this result gives back exactly Eq.(47l for a^ = /(O), which shows the 



It is possible to write Eq. (k34| in a simpler form at the price of a small approximation, by making use of the Stirling 
expansion of the Euler function Inr(z) ~ (z - 1/2) Inz - z + (1/2) ln(27r). Writing r^{ix) = \T{ix)\'^ {r{ix)/T{-ix)) 
where [22] |r(ia::)p = tt/[x sinh(7ra;)], with r{ix)/r{—ix) = — r(l + ix)/r{l — ix) and making use of the Stirling 
expansion for r(l ± ix), we find: 



/(^) 



sinh(7ra;) 



sm 



„ , , 9x arctanx 
6 + 2 - ln(l + x^) I X 



(65) 



with the notation: 



b :=: \n{m^ / m^^) — ln(ln(?n-|-/?7i^)) 



(66) 



This result is an excellent approximation for Eq.(64), except in the vicinity of a; = because it gives /(O) = b + 1 
instead of /(O) = b + 2C, which is still a very good approximation since 2C ~ 1.154. This small deficiency comes 
from the Stirling expansion in the vicinity of z = 1. We may correct for it by replacing 2 by 1 + 2C. This leads to: 



/(^) 



sinh(7rx) 



6 + 1 + 2C - ln(l 



arctan x 



(67) 



which is almost indistinguishable from Eq.(64), for any sizeable val ue of b. Thi s ca n be seen from Fig. p^ where, for 
m^/ra^ — 10^ corresponding to & ~ 3.08, we have plotted both Eq.(64) and Eq.(67). They are undistinguishable on 



the figure. We have also plotted the exact solution of the first order approximation Eq.(31 ), for small m^/m^, of the 
exact starting equation Eq.([3|. The agreement is very good. The only sizeable difference is again in the vicinity of 
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FIG. 4: (color online). Integral of Eq.(67l as a function of 6 = ln{m^/m_i) — ln(ln(m^/m^)) compared to the exact result 7r/2 
from the sum rule Eq.([5|. 



X = because the asymptotic form Eq.(47) a^ — b + 2C is, for this value of m^/m^ ~ 10^, still slightly different 
from the result obtained from Eq.(31 1 (or equivalently Eq.(37)) and Eq.(|3|, as it is seen from Fig. IT] This could be 
completely corrected by replacing, in Eq.(67l, & + 2C by the exact 03. This would provide an excellent approximation 
to the exact solution of the original Eq. ( SFT Finally we have also checked that this expression Eq. ( [67| satisfies the 
sum rule Eq. ([5]) with a remarquable precision as it can be seen from Fig. HI 

This explicit result Eq.(67 1 shows that f{x) decreases very rapidly, with an oscillatory behaviour in an exponentially 
decreasing envelope. This function is almost zero after the first zero, except for a small oscillation just after this zero. 
For large 6, this zero is small and is given essentially by xq = 7r/(6 + 2C) = tt/Ss. Hence we find explicitely a 
result coherent with our starting hypothesis, namely that f{x) gets very narrow for large 03. Actually the behaviour 
resulting from Eq.(64) for extremely large x is not the correct one because we have made approximations for the low k 
behaviour of f{k). As we mentioned above this behaviour is rather f{x) ~ l/x^'^^''^ without any oscillation. However 
since f{x) is anyway extremely small in this regime, this is completely unimportant in practice. We note that all our 
above analytical findings are in excellent agreement with the very recent numerical work of Iskin [17]. A final point, 
also in agreement with the work of Iskin [T7] is that the first zero disappears (by coalescence with the second one) 
when b decreases. This feature is fairly natural since, when m^ goes to infinity, the solution is the lorentzian found 
in section HI which has no zero. Indeed the argument of the sinus in Eq.(67) is an increasing function of x for small 
X, but decreases for large x. The first zero disappears when this argument reaches its maximum as a function of x 
and is equal to tt at this maximum. This occurs for Xq = 2.33 and b = 1.55, corresponding to m^/m^ — 11.6. This 
feature is again in very nice agreement with the full numerical solution [17 of Eq.(|3]), although it is not very good 
quantitatively since the zero is found to disappear for m-^/m^ ~ 7 corresponding to b =^ 1.28. This small quantitative 
disagreement is expected since this mass ratio does not correspond to the full asymptotic regime of very small m^. 

We note finally that, if we are not interested in the disappearance of the zeros and rather concentrate on the regime 
where 03 is quite large, f{x) can be further simplified into a very simple expression. Indeed we may assume that x 
is fairly small since otherwise, from Eq.(67l, f{x) is extremely small, compared to its a; = value, and its specific 
expression becomes unimportant. 
03 = & + 2C, Eq.([67| reduces to: 



In this case we have ln(l + x ) ~ and arctana:/x ~ 1. Taking into account 



/(^) 



sinh(7ra;^ 



sin(a3a;) 



(68) 



First this expression satisfies /(O) — 03. Then with this very simple form the integral in the sum rule Eq.(l5]) can be 
calculated analytically. One finds {tt/2) tanh(a3/2) which, for large 03, is 7r/2 as it should be. Furthermore when the 
right-hand side of Eq.Q is calculated analytically one finds, to dominant order, the result 03 e~°'^. This is in perfect 
agreement with the left-hand side of Eq.(K3|) since /(O) = 03 and, to dominant order, R/2 = r = m^/ra^ ~ e~°'^ from 
Eq.(47|. Hence, since a^ is large, we have an almost complete cancellation of the two terms in the right-hand side of 
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Eq. (p| whereas a simple order of magnitude evaluation would rather give a comparatively much larger result of order 
(1/53)^. It is thus quite remarkable that the non trivial property Eq.(k3| is satisfied by the simple expression Eq.(68|. 



V. CONCLUSION 

In this paper we have considered the problem of the scattering length 03 for a fermion of mass m^ colliding 
with a dimer, formed from a fermion identical to the incident one and another different fermion of mass ?7i^. The 
only scattering parameter in this problem is the scattering length a between the different fermions. We have been 
interested in the way in which this scattering length depends on the mass ratio m^/m-^ between the different fermions. 
The experimental investigation of these kind of systems is under current development in several laboratories in the 
world. While the answer to this problem is obtained from a generalization of an integral equation first explored by 
Skorniakov and Ter-Martirosian [14j (STM), which has been fully investigated numerically, we have been interested 
in an analytical investigation of this problem and this equation in order to gain additional insight and open possible 
roads for more complex problems. 

We have been able to perform this investigation in the limiting cases where the mass ratio is very large or very small. 
When the mass m^ of the lonely fermion is very large, the situation is very similar to the one where it is infinite and 
where the corresponding atom behave as a fixed scattering center. In this case we have found the analytical solution 
of the STM integral equation. The corresponding scattering length is merely a. The existence of another fermion in 
a bound state (i.e. the existence of the dimer) is unimportant and physically everything happens as if the fermion 
is merely scattering on the fixed center. It is noteworthy that the equal mass case m^ — m^, where a^/a ~ 1.18, is 
fairly near this limiting case. 

When the mass tti^ is very small, the atom-dimer scattering length a^ goes to infinity. We have found that in 
this situation the width of the unknown function in the generalized STM equation becomes very narrow around the 
origin. We have used this feature to expand properly the STM equation in this case. We have in this way obtained 
the explicit asymptotic behaviour of 03 in this large m-^/m^ regime. We have found: 

— ~ ln(m^/TO^) - ln(ln(TO^/TO^)) + 2C (69) 

where C is the Euler constant. In addition we have found a parametric representation of 03 in terms of m^ /m^ 
which is even more accurate at lower m^/mi, and gives an excellent agreement with the exact numerical results over 
the whole range m^/m\^ > 7. At the same time we have found an excellent approximate analytical solution for the 
unknown function in the STM equation in this same regime. In particular this solution displays the same appearance 
of zeros and oscillations as it is found in the numerical results of the exact equation for m-f/mj^ > 7. This is in contrast 
with the shape of this solution for lower m-^/m^ which is merely found to decrease monotously to zero. 

Appendix A: Interpretation of the integral relation EqJ5] 

In this Appendix, we show in details the connection between the integral relation Eq. ([5]) and the fact that two spin 
up particles can not be at the same position. In order to do so, we need to make the connection between the function 
o-siP) of Eq.(l2]) and the wave function /(r) studied by Petrov ^8j. We keep his notations, so we emphasize that in 



this Appendix, contrary to the rest of the article, / is a wavefunction defined in Eq.( Al ) below. 

We first recall the notations and results of il8j. The system we study is made of two spin up particles (positions 
Xi and X2) and a single spin down particle (position X3). We want to describe the scattering between a spin up 
particle and a dimer in the limite of zero incident energy. The total energy E is therefore equal to the binding energy 
E = — l/(2/Lta^) < 0. In the limit where the spin down particle of position X3 is close to the spin up particle of 
position X2, the wave function Vt is given by 

*^M(_L_^) (Al) 

where rj^ = 2 ^Jm-^ (X3 — X2). In the limit X2 — > X3, the vector r is given by r = y'm^(x2 — Xi)/sin0 (therefore r 
has the dimension of a length multiplied by the square root of a mass, or equivalently, of the inverse square root of 
an energy, since we take h = 1). The angle 9 is given by tan^ = y/1 + 2m^/m^. In Ref.fTQ, the following linear 
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homogeneous equation for /(r) was derived 

4^ / dr' [CEi^n 



r'2 -2r-r' cos(26l))/(r') 
+ G^(|r-r'|)(/(r)-/(r')) 



(A2) 



with GfidXI) = J jl^ ^* ^ "^ TTeI+ioT ■ ^^ '^^'^ Fourier transform Eq.(|A2|) and we find 



d3q' 



/(q') 



(27r)3 [|£;| sin^ + (q^ + q'2 - 2 cos 61 q • q') sin6l] 
1 



47r 



^/\El^-^\ /(q) = 



Next, we define the function g sucli that 



4-7r 
/(q) ^ (2^)35(q)--7 5(7q) 



(A3) 



(A4) 



where 7 is a constant which is unknown a priori. Inserting Eq.(|A4| in Eq.(A3), we find that we must take 7 



2 a ^/m^/^/l + 2 r a-^/m i in order to recover Eq.|2|. We thus have ^(q) = a3(q). 

Integrating Eq.(A4) we can now express the integral relation Eq.(l5]) for 03 in terms of the function /(r) 



d^q 



/(q) 



(27r)3- 
-- fir = 0) 



(27r)3^ 
2 '■°° 



75(7q) 



dqa3{\q\) 



(A5) 



However, r = corresponds to the situation where the two spin up particles are at the same position (xi — X2) 
and therefore the wave-function ^, hence /, must vanish due to the Pauli exclusion principle : /(r = 0) = 0. From 
Eq.(A5), we see that we must have 



dqasilql 



■K 

2 



(A6) 



which is precisely Eq.([5]). As a consequence, we see that the integral relation Eq.([5]) is simply equivalent to the 
statement that two spin up particles can not be at the same position in space. This also means that it is not expected 
to be correct for three identical bosons, for instance. 
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